Structure formation in the DGP cosmological model 
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The DGP brane-world model provides an alternative to the standard LCDM cosmology, in which 
the late universe accelerates due to a modification of gravity rather than vacuum energy. The 
cosmological constant A in LCDM is replaced by a single parameter, the crossover scale r c , in DGP. 
The Supernova redshift observations can be fitted by both models, with A ~ Hq and r c ~ Hq ■ This 
degeneracy is broken by structure formation, which is suppressed in different ways in the two models. 
There is some confusion in the literature about how the standard linear growth factor is modified 
in DGP. While the luminosity distance can be computed purely from the modified 4-dimensional 
Friedman equation, the evolution of density perturbations requires an analysis of the 5-dimensional 
gravitational field. We show that if the 5-dimensional effects are inappropriately neglected, then 
the 4-dimensional Bianchi identities are violated and the computed growth factor is incorrect. By 
using the 5-dimensional equations, we derive the correct growth factor. 



I. INTRODUCTION 

The acceleration of the late-time universe, as im- 
plied by observations of Supernovae redshifts, cosmic 
microwave background anisotropies and the large-scale 
structure, poses one of the deepest theoretical problems 
facing cosmology. Within the framework of general rela- 
tivity, the acceleration must originate from a dark energy 
field with effectively negative pressure, such as vacuum 
energy or a slow-rolling scalar field ( "quintessence" ) . So 
far, none of the available models has a natural expla- 
nation. For example, in the simplest option of vacuum 
energy, leading to the "standard" LCDM model, the in- 
credibly small, 
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and incredibly fine-tuned, 
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value of the cosmological constant cannot be explained 
by current particle physics. 

An alternative to dark energy plus general relativity 
is provided by models where the acceleration is due to 
modifications of gravity on very large scales, r > Hq 1 . 
One of the simplest covariant models is based on the 
Dvali-Gabadadze-Porrati (DGP) brane-world model Q, 
in which gravity leaks off the 4-dimensional Minkowski 
brane into the 5-dimensional "bulk" Minkowski space- 
time at large scales. At small scales, gravity is effectively 
bound to the brane and 4-dimcnsional Newtonian dy- 
namics is recovered to a good approximation. The tran- 
sition from 4- to 5-dimensional behaviour is governed by 
a crossover scale r c ; the weak-field gravitational potential 
behaves as 



r 1 for r < r c 
r~ 2 for r > r c 



(3) 



The DGP model was generalized by Deffayet to a 
Friedman-Robertson- Walker brane in a Minkowski bulk 



; the gravity leakage at late times initiates acceleration 
- not due to any negative pressure field, but due to the 
weakening of gravity on the brane. The energy conserva- 
tion equation remains the same as in general relativity, 
but the Friedman equation is modified: 



SH(p + p) =0, 
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It is important to stress that the modification to 
the Friedman equation is derived from a covariant 5- 
dimensional action and junction conditions across the 
brane Q. 

The modified Friedman equation (J5J shows that at late 
times in a CDM universe, with p oc a -3 — > 0, we have 



H 



(6) 



Since Hq > Hoo, in order to achieve acceleration at late 
times, we require r c > H^ 1 , an d this is confirmed by 
fitting SN observations Q. 

Like the LCDM model, the DGP model is simple, with 
a single parameter r c to control the late-time accelera- 
tion. The DGP model does not provide a natural so- 
lution to the late-acceleration problem; similarly to the 
LCDM model, where A must be fine-tuned, the DGP 
parameter r c must be fine-tuned to match observation. 
Furthermore, it has been recognized that the DGP model 
suffers from potentially serious theoretical problems, such 
as the existence of a ghost in de Sitter solutions of Eq. (J5J 
and the strong coupling problem ,4]. Nevertheless, like 
LCDM, the DGP is a simple covariant model that makes 
testable predictions, in particular the prediction of late- 
time acceleration. (It is worth noting that the original 
DGP model with a Minkowski brane was not introduced 
to explain acceleration - the generalization to a Friedman 
brane was subsequently found to be self-accelerating.) 

LCDM and DGP can both account for the SN observa- 
tions, with A ~ Hq and r c ~ Hq 1 respectively. This de- 
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generacy may be broken by observations based on struc- 
ture formation, since the two models predict different 
suppression of density perturbations. For LCDM, the 
analysis of density perturbations is well understood. For 
DGP it is much more subtle and complicated; although 
matter is confined to the 4-dimensional brane, gravity 
is fundamentally 5-dimensional, and the bulk gravita- 
tional field responds to and backreacts on density per- 
turbations. The background cosmological dynamics, and 
hence the predicted SN redshifts, can be treated purely 
via the modified 4-dimensional Friedman equation. By 
contrast, the evolution of density perturbations requires 
an analysis based on the 5-dimcnsional nature of gravity. 



II. PERTURBATION EQUATIONS ON THE 
BRANE 



Effective covariant equations on the brane, for an ar- 
bitrary brane metric and matter distribution, can be de- 
rived by projecting the 5-dimensional Einstein equations 
(Gab = 0) an d applying the Israel junction conditions 
with reflection symmetry at the brane 0, |(j 



= (167rGr c ) 2 n M „ - E, 
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and E^v is the trace-free projection of the 5D Weyl ten- 
sor. The energy-momentum tensor satisfies 



= 0, 



(10) 



as in general relativity. The 4D Bianchi identity imposes 
the constraint 



y^E^ = (167rGr c ) 2 VTV 



(11) 



In the background spacetime, E^ v — and Eq. Q 
implies the modified Friedman equation JSJ. The per- 
turbed FRW brane has a nonzero £^„, which encodes 
the effects on the brane of the bulk gravitational field. 
The perturbed 5D field equations are needed to deter- 
mine the evolution of SE^. We can parametrize the 
(gauge-invariant) scalar perturbations of E^ v as an ef- 
fective fluid, with density perturbation SpE, isotropic 
pressure perturbation j^SpE, anisotropic stress pertur- 
bation Site and energy flux perturbation SqE (see the 
appendix). Scalar metric perturbations are given in lon- 
gitudinal gauge by 



and the perturbed energy-momentum tensor for matter 
is given by 



-Sp 
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The Bianchi identity pi(l gives evolution equa- 
tions HA2|1 and (|A3|) for SpE and SqE- The perturbed 
4D field equations Q imply a modified Poisson equation 
(IA4() and constraint equation Ijltifl . The 4D perturbation 
equations are not closed, due to the presence of Weyl 
anisotropic stress Site- The 5D perturbation equations 
are needed to determine this stress. This is an extremely 
complicated problem in full generality, and especially on 
large scales where the 5D effects become strong. For 
the structure formation problem, we focus on subhorizon 
scales, with k/a -C H^ 1 and k/a <C r c , where k is the co- 
moving wavenumber. In addition, we use the quasi-static 
approximation which is relevant for structure formation 
and neglect time-derivatives terms relative to gradient 
terms. Then Eq. I|A2[) implies 



Sq E = 0. 



(14) 



With SqE — 0, the key perturbation equations reduce 



to 



K* = ^(J^)(pA-£-), (15) 
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where 
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\-2Hr c r c (2H 2 + H)-H 2Hr c -V 

(19) 

and the (gauge-invariant) comoving density perturbation 
is 



pA = Sp — SHaSq . 



(20) 



The equations (|15fl - i|18[l are not closed since SE^ is 
not determined by purely 4D equations. The ad hoc as- 
sumption SE^ — 0, i.e., 

Sp E =Stt e = 0, (21) 

has been effectively adopted in some papers Q. However, 
it is manifest that this assumption is not consistent with 
the Bianchi identity Eq. i|17|) . Thus the analysis based 
on this assumption must be revisited. 

Within our approximations, the Bianchi identity sug- 
gests that there exists a family of solutions characterized 
by 



ds 2 



-(1 + 2^)dt 2 + a 2 (l + 2$)dx 2 



(12) 



Sp E = C{Hr c )k 2 6ir E ■ 



(22) 
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where C(Hr c ) is some function of Hr c . The correct 
C(Hr c ) can be obtained only when we properly consider 
the 5D equations of motion. 



III. THE 5-DIMENSIONAL EQUATIONS 

The 5D background Minkowski metric is given by 



dsf 5) = -n(y, t) 2 dt 2 + b(y, tfdx 2 + dy 2 
where the brane is at y = and 

b = a(l + Hy), n = 1 




H y. 



(23) 



(24) 



Because the background bulk is just a Minkowski space- 
time, we can easily find a Mukohyama master vari- 
able tt(t,x,y) from which all the solutions for the 
metric perturbations can be constructed (see the ap- 
pendix). In the quasi-static approximation relevant for 
structure formation, we can neglect Q relative to f2', and 
the 5D wave equation for f2 becomes 



4 + «nt 

b n 



k 2 
b 2 



n = o, 



(25) 



where a dot denotes a derivative with respect to time t 
and a prime a derivative with respect to y. In the quasi- 
static approximation, the solution is 



n = ci 



(1 + H V y k / aH + c 2 (l + Hy) k / aH 



(26) 



Regularity of the bulk perturbations as y — ► oo requires 
that C2 = 0. 

Perturbations of the Weyl fluid SE^ V can be expressed 
in terms of Q, 
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y=o 



By Eq. (|A14|) with C2 = 0, we can neglect the last term 
on the right of Eq. . It follows that 



SpE = 2fc She , 



(29) 



and we have determined the function C in Eq. I|22ll as 
C = 2. 



using Eqs. ||aT|) . (|A8|) . (jATTI) . {HZ) and (|A14|) . the so- 
lutions for the brane metric perturbations are 



5* = 4 ^( 1 "P)^ 
5* = -4^(l + ^)pA 



where 



/3 = 1 - 2r c H 1 



H 
3H 



2 I ' 



(30) 
(31) 

(32) 



This agrees with the results obtained by Lue, Scocci- 
marro and Starkman 9]. They find spherical symmetric 
solutions by closing the 4D equations using an anzatz for 
the metric and checked in retrospect that the obtained 
solutions satisfy the regularity in the bulk. Here we have 
shown that the solutions H3Ufl and l|31|l are uniquely de- 
termined by the regularity condition in the bulk within 
our approximations. 

We can check that our solutions are consistent with the 
Bianchi identity. The Bianchi identity Eq. (|17|) gives non- 
trivial equations for 8pE and Site with the condition l|29() 
and they can be solved in terms of pA. Then substituting 
the solutions into Eqs. and (|lf>|> . we arrive at the 
solutions for $ and * as in Eqs. (J20J) and (J3TJ. 

The modified Poisson equation (|30l) shows the suppres- 
sion of growth due to gravity leakage. The rate of growth 
is determined by A, and for CDM, 



k 2 

2HA = — . 

a 2 



which leads to 



A 



A 



2HA = 4ttG I 1 + — 



P A 



(33) 



(34) 



Thus the growth rate receives an additional modification 
from the time variation of Newton's constant through (3. 
This effect can be described by a linearized scalar-tensor 
gravity with Brans-Dicke parameter Q 



0J = -(f3-l), 



(35) 



where the gravitational scalar corresponds to the bending 
of the brane Eq. I|A9|I . 

In Fig. 1, we show the linear growth factor A/a for the 
DGP model, and compare it with LCDM, with the incor- 
rect DGP result (in which the inconsistent assumption 
SE^ = is effectively adopted), and with the general 
relativity dark energy model whose background evolu- 
tion matches that of the DGP model (see Fig. 2). 



V. CONCLUSION 

IV. THE GROWTH FACTOR 

We have shown by a careful analysis of the 5D per- 
In the appendix, we find the solution for l|A15|) by turbation equations that subhorizon density perturba- 
imposing the perturbed junction condition (|A13(1 . Then tions in the DGP cosmological models evolve according 
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FIG. 1: The growth history g(a) — A(a)/a is shown for 
LCDM (long dashed) and DGP (solid, thick). The growth 
history for a dark energy model (short dashed) with the iden- 
tical expansion histories with DGP is also shown (see Fig. 2) 
[TrT |. Due to the time variation of Newton's constant through 
/3 in Eq. (JHU, the growth factor g(a) receives an additional 
suppression compared with the dark energy model. DGP-4D 
(solid, thin) shows the incorrect result in which the incon- 
sistent assumption SE^v — is adopted. We set the desity 
parameter for matter today as f2 m o = 0.3. 



to Eq. I|34|). that the modified Poisson equation is given 
by Eq. (|30|l . and that the gravitational potential is given 
by Eq. I|31(l . These expressions obey the constraint im- 
posed by the 4D Bianchi identity, and they correct those 
expressions in the literature |jj which violate the Bianchi 
identity. Our results confirm the results derived via an- 
other approach by Lue, Scoccimarro and Starkman 0. 
Although they also used the 5D equations to check con- 
sistency of their ansatz, they were not able to show, as 
we have done here, that their results follow uniquely from 
the 5D equations. 

The correct equations for subhorizon density perturba- 
tions are crucial for meaningful tests of DGP predictions 
against structure formation observations. And such tests 
are essential for breaking the degeneracy with LCDM 
that arises with SN redshift observations 0, El- The 
distance-based SN observations draw only upon the back- 
ground 4D Friedman equation JJJJ in DGP models, and 
therefore there are quintessence models in general rela- 
tivity that can produce precisely the same SN redshifts 
as DGP (see Fig. 2). By contrast, structure formation 
observations require the 5D perturbations in DGP, and 
one cannot find equivalent general relativity models. Pre- 
vious analyses of structure formation that are based on 
the incorrect equations for density perturbations |7| will 
need to be re- visited. 

The linear growth factor, as shown in Fig.l, is the ba- 
sis for tests of DGP against structure formation observa- 
tions. But one also needs the metric perturbations. It 
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FIG. 2: The expansion history r(z) — J* dzH(z) 1 is shown 
for LCDM (dashed) and DGP (solid) where 1 + z = o" 1 . 
Expansion history for DGP can be mimicked by a dark energy 
model with the equation of state w(a) — wq + w a (l — a) with 
wo = -0.78 and w a = 0.32 []]J. 



follows from Eqs. and that 
k 2 

($_$) = SirGpA, (36) 

a z 

k 2 , , 8ttG a 

- $ + * = ——pAocSwE, 37 

a z 3p 

where the proportionality follows from Eq. Equa- 
tion (|37|l differs from general relativity, where the ab- 
sence of matter anisotropic stress implies ^ = — $. In 
the DGP model, the 5D gravitational field produces an 
anisotropic stress on the 4D universe. Equation (|36|l is 
the same as in general relativity - and this equation ba- 
sically determines the integrated Sachs- Wolfe and weak 
lensing effects. This means that the formulas will take 
the same form as in general relativity, but the outcome 
will be different, because A will differ from the general 
relativity case. 

We have also emphasized that our results (and those 
of Ref. 0) apply on subhorizon scales. On superhori- 
zon scales, where the 5D effects are strongest, the prob- 
lem has yet to be solved. Thus a computation of CMB 
anisotropies will not be reliable on the largest scales. 

Furthermore, we have restricted our discussion to the 
linear regime. In the nonlinear regime, the DGP theory 
approaches general relativity ffl| . A full computation of 
weak lensing will require a careful matching of the linear 
DGP to the nonlinear general relativity regimes. 
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APPENDIX A: PERTURBATION EQUATIONS 

We can parametrize the scalar perturbations of as 
an effective fluid, 



a 1 dq E 



aSqs,i 



3 •■ S *Ej 



(Al) 



where Snf- — 5ir^ 
Eq. l(TT|) imposes 



• ■• 1 1 -\""E k u ij- The 4D Bianchi identity 



\5ir 



8p E + AH8p E -a 1 k 2 5q E = 
Sq E + AH5q E + " 



(A2) 



1 2 



,2 H 
= —a — r r — 
3 H 

1 



P A 



Sp E 



r c (2H + H/H) - 1 



2Hr c - 1 2 J ffr„ 
k 5tt e 



(A3) 



Then the induced metric perturbations on the brane are 
given by 



q> = A- I^ + hU, $ = (A10) 



from a gauge transformation. We can express ^ and $ 
in terms of A and R as 



$ = 



* = 




(All) 



1-r, 



(l-Hr c )A + r c H+ „ H2 R ] , 



H 



(A12) 



Combining the (i, f) component and the (0, i) component 
of the perturbed 4D field equations Eq. Q , we get the 
modified Poisson equation 



— $ = AttG 1 



\2Hr 



pA 



SpE — SHaSqE 



(A4) 

where pA = 5p — 3Ha5q. The traceless part of com- 
ponent gives Eq. (JTBJ. (We follow the standard assump- 
tion that the CDM anisotropic stress vanishes Stt = 0.) 

In order to solve the perturbations in this background, 
it is convenient to use 5D Longitudinal gauge, given by 



ds 2 = (l + 2A yy )dy 2 + nA y dydt - n 2 (1 + 2 A 2 )dt 2 



-b 2 (l + 2K)dx 2 . 



(A5) 



In Minkowski spacetimc, if the metric perturbations are 
derived from a "master variable", fl, the perturbed 5D 
Einstein equations yield a single wave equation governing 
the evolution of the master variable Q in the bulk: 



1 

nb 3 



n 



V«)'-lF a = ' (A6) 



Within our approximation, we can neglect time- 
derivative terms and the metric perturbations are written 
in terms of Q as 



.4 
R 



1 

66 



b 2 



(A7) 
(A8) 



An important feature of the DGP model is the brane- 
bcnding mode; in 5D longitudinal gauge, the location of 
the brane is perturbed and given by |12( 



V = £ = -r c ($ + ¥). 



(A9) 



The perturbed junction condition is 



k 2 k 2 P 1 k 2 



2a 2 r r 



4r a 3 l v =o- 



(A13) 



This junction condition gives the modified Poisson equa- 
tion where the last two terms describe corrections from 
the 5D gravity. 

Note that fl' appears in the expressions for SpE and 
8-ke- Unless we solve the 5D equations and impose proper 
boundary conditions, fl and f2' on the brane are inde- 
pendent. This means that we can impose arbitrary con- 
ditions for the relation between SpE and Site- Indeed 
by choosing appropriate c 2 in Eq. I|26l) . we can have any 
C(Hr c ) in Eq. (|22l) . However, as described in the main 
text, the regularity condition impose c 2 = 0. Then we 
get 



k 2 

HQ' < —zQ, 



(A14) 



for k I aH 3> 1 and we can neglect the terms proportional 
to il'. However, the Poisson equation Eq. (|A13|I still re- 
ceives a correction from the brane ben ding mode £ and we 
get a linearized scalar-tensor gravity 13]. The junction 
condition Ea. (|A13|l gives 



k 4 Q _ 1 - 2Hr c (l + H/2H 2 ) 
2^~ ~ i _ 2Hr c (l + H /3H 2 ) 



nip A. 



(A15) 



Then from Eqs. IpTF)). iJXSJl, (fXTTj) and JXTSJ, the so- 
lutions for metric perturbations can be obtained as 
Eqs. H3U|) an d (|31|l . We can also derive the solutions 
for SpE and Site from Eqs. I|27|l and (|28|) and show that 
they are consistent with the 4D Bianchi identity. 
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